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Abstract 

We analyze general structure of AA-fold supersymmetry which provides a systematic framework 

to construct weakly quasi-solvable quantum mechanical systems. Main ingredients of our analysis 

are dimensional analysis and introduction of an equivalent class of linear differential operators 

associated with A/"-fold supersymmetry for each M. To illustrate how they work, we construct the 

,4^ | most general form of TV-fold super symmetric systems for J\f = 2, 3, and 4. 
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I. INTRODUCTION 

A/"- fold supersymmetry (SUSY) [1-3] is one of the most powerful frameworks for construct- 
ing a one-dimensional quantum mechanical (QM) system which admits analytic solutions 
in closed form in a certain sense. This is due to the fact that A/-fold SUSY is essentially 
equivalent to weak quasi-solvability which is until now the least restrictive concept about the 
availability of solutions in closed form. The latter crucial fact was first proved in a general 
fashion in Ref. [2] and was later complemented slightly in Ref. [4]. For a review, see, e.g., 
Ref. [5]. 

In general, construction of an A^-fold SUSY system get more difficult as the number 
M e N increases since we must solve coupled nonlinear differential equations for A^ unknown 
functions (see the next section). To bypass the latter difficulty, a systematic algorithm for 
constructing an TV-fold SUSY system based on quasi-solvability (in the strong sense) was 
proposed in Ref. [6]. A key ingredient of the algorithm is to choose first an A/"-dimensional 
linear space of specific functions such that it can be preserved by a second-order linear 
differential operator. It has been proved to be quite efficient and so far four inequivalent 
types of A/"-fold SUSY, namely, type A [4, 7], type B [8], type C [6], and type X 2 [9] are 
successfully constructed with the algorithm. We note that almost all the models having 
essentially the same symmetry as A/"-fold SUSY but called with other terminologies in the 
literature, such as Poschl-Teller and Lame potentials, are actually particular cases of type 
A TV-fold SUSY. 

It is evident, however, that the algorithm is helpless to construct a weakly quasi-solvable 
system which only admits a finite-dimensional invariant subspace determined by another 
differential equation. The framework of AA-fold SUSY covers such systems and thus provides 
a more general formalism than higher- derivative generalizations of Darboux transformation 
such as the Crum's method [10] which relies on a set of exact eigenfunctions of a regular 
Sturm-Liouville system. To construct a weakly quasi-solvable system, we must in general 
treat directly the aforementioned coupled nonlinear differential equations. For the simplest 
case of M = 2, the general result was already studied and reported in Refs. [2, 11, 12]. On 
the other hand, for the cases of J\f > 2 until now on there is, at the best of our knowledge, 
only one paper [13] which studied the M = 3 case. This fact would reflect the difficulty and 
complexity of the problems for larger A/". 

In this work, we investigate general strucuture of A/"- fold SUSY systems to extract relevant 
clues to construct them which have in particular weak quasi-solvability. For this purpose, we 
first employ dimensional analysis which is a well-known powerful tool in general physics. It 
turns out that it is also quite efficient in acquiring deeper understandings of A/"-fold SUSY. 
We then introduce equivalent classes of linear differential operators associated with AA-fold 
SUSY. We find that it enables us to deal with operator equalities appeared in A/"-fold SUSY 
more systematically and transparently. 

We organize the paper as follows. In the next section, we first briefly review the in- 
gredients of AA-fold SUSY. Then, we introduce two key concepts for analyzing its general 
structure, namely, dimensional analysis and equivalent classes of linear differential operators 
associated with AA-fold SUSY. In Sections III— IV, we apply the general arguments to obtain 
general form of A/"-fold SUSY for M = 2, 3, and 4, respectively. We show how dimensional 
analysis enables us to reduce the complexity of the problems on solving the conditions for 
A/"- fold SUSY and on finding integral constants of the systems. In the last section, we 
summarize the paper and provide comments on the future issues. 



II. GENERAL CONSIDERATION 

To begin with, we shall briefly review ingredients of jV-fold SUSY as preliminaries. For 
details, see the review [5]. An A^-fold SUSY QM system in one-dimension is composed of a 
pair of Hamiltonians H ± and a pair of A/th-order linear differential operators Pj^ 



w+^>> *-£,+%<**>& 



H± = -oh + V± (1)< P » = hi + £ ^Wfr. ^=TO) T . (2.1) 



fc=0 



where the superscript T denotes the transposition of a linear operator [3], which satisfy the 
intertwining relation 



fc=0 



P^H-^H-P x = J2l^ k =0, (2.2) 



and its transposed relation H~P^- — P^-H + = 0. The operators P^- are actually components 
of A/"-fold supercharges. 

One of the most significant consequences of the intertwining relation (2.2) is weak quasi- 
solvability [2, 4]. That is, each A/"-fold SUSY Hamiltonian H^ preserves the linear space 
ker Pp. 

/^kerP^ckerP^. (2.3) 

If the differential equation P^<p = and/or P^j-(f> = admits a number of analytic solutions 
in closed form, P~ and/or H + is not only weakly quasi-solvable but also quasi- solvable in the 
strong sense. But in general, an A/"-fold SUSY Hamiltonian is merely weakly quasi-solvable 
and does not admit any analytic local solutions. We also note that ker P^ is not necessarily 
a subspace of the linear space, which is usually the Hilbert space L 2 (S) (Set), in which 
the operator H ± acts. 

Another peculiar feature of an A^-fold SUSY system is that the product PjfP^f which 
arises as a component of the anti-commutator of A/"-fold supercharges is an AAth-degree 
polynomial in the Hamiltonian H ± [2, 3] and thus has the following form: 



r N r M — z 



Af-l 

M-k-l 



(H ± + Co) N +Y. C ^ H± + C » 



fe=i 



(2.4) 



where Ck {k = 0, . . . , M — 1) are all constants. The M zeros of the polynomial in the 
r.h.s. of (2.4) correspond to the spectrum of H ± in the space kerP^. Hence, they are 
actually a part of the eigenvalues of H ± if kerP^ C L 2 (S). In the latter case, we can 
calculate the part of the eigenvalues algebraically from (2.4) even though the corresponding 
eigenf unctions cannot be obtained in closed form. We note that the M energy spectra would 
exhibit the characteristic features such as the disappearance of nonperturbative corrections 
due to the generalized non-renormalization theorem in A/"-fold SUSY, which were reported 
on the realistic physical systems such as the asymmetric double-well potential in [14] and 
the periodic and the symmetric triple- well potentials in [15]. 

One of the most difficult problems on A/"-fold SUSY is to analyze the condition (2.2) 
for A/"-fold SUSY. It is composed of coupled nonlinear differential equations for the so far 



undetermined functions w k (k — 0, . . . , M — 1). As is easily expected, its complexity gets 
terrible as the integer M increases. Hence, it is quite difficult to solve directly the condition 
(2.2) for larger N ' . In the subsequent two sections, we shall investigate general aspects of 
the complicated structure of A/"-fold SUSY systems which would provide us a clearer view 
on them. 



A. Dimensional Analysis 

Dimensional analysis is one of the powerful methods to make physical consideration and 
in particular to estimate a physical quantity under consideration without solving equations 
directly, cf., any textbook on general physics. In this section, we shall see that it supplies 
us with a valuable guiding principle in solving the condition for A/"-fold SUSY. 

To make dimensional analysis on our system (2.1), we first note that we have implicitly 
employed the unit system where the ratio of square action (the Planck constant) to mass 
H 2 /m is dimensionless. The only relevant physical dimension is then the length, denoted by 
[L], which is carried by the physical position variable q. It is easy to see from (2.1) and (2.4) 
that the physical dimensions of V ± , C&, and w k in terms of the length are given by 

V ± [L~ 2 ], C k [L~ 2 ^] (k = 0,...,Af-l), 

w \mm) [L k^- m] (A; = 0, . . . , AT - 1, m = 0, 1, 2, . . .), ^"^ 

where w k (q) is the mth derivative of w k (q) with respect to q. 

Dimensional analysis relies on the obvious fact that all the terms which appear in a 
single formula under consideration must have the same physical dimension. For instance, a 
potential has the physical dimension [L~ 2 ] and thus must be expressed as a sum of terms all 
of which have the same physical dimension [L -2 ] . Hence, if we only consider a polynomial 
of Ck and w k (m — 0, 1, 2, . . .), it must have the following form 

V = a w [ f_ 2 + axwjJO; + a 2 (4^-i) 2 " Co [L~ 2 ], (2.6) 

where a k (k = 0, 1,2) are all dimensionless parameters. In fact, we can see that a pair of 
JV-fold SUSY potentials V ± satisfying (2.2) does have the form (2.6). The l.h.s. of (2.2) 
is a linear differential operator of at most A/th order, as is indicated in (2.2), and it is 
evident that the identity (2.2) holds if and only if all the coefficients I k of d k = d k /dq k 
(k — 0, . . . , Af) vanish. The latter requirement for the coefficients of d and d reads as 

^ 1 =^'i-(^ + -^)=0[L- 2 ], (2.7a) 

2$!i = "#5 + 2^jf 2 + 2AA/-' - 2<1 X (U+ -V~) = [L- 3 ]. (2.7b) 

The set of conditions (2.7) can be easily solved as 



V ± - —*>&* + "^ ± -o) -K + ^7«^) 2 - Co [L- 2 ], (2.1 



M w ^ 2+ \ 2M ± 2) WN ~ l + 2M [ 

which indeed has the form of (2.6). The formula (2.8) provides a general expression for a 
pair of A/"-fold SUSY potentials V^ for an arbitrary Af 6 N. One of its characteristic features 



is that they are expressible solely in terms of the two functions Wj^^ and wjv_ 2 irrespective 

of what additional conditions w k (k = 0, . . . , TV — 1) should satisfy. 

The remaining conditions for TV-fold SUSY coming from the coefficients of d k for k = 
0, . . . ,TV — 2 in (2.2) are in general algebraic equations consisting of w k (m = 0, 1,2, . . .) 

after the substitution of (2.8) into (2.2). Dimensional analysis tells us that the operator in 
the l.h.s. of (2.2) has the physical dimension [L ~ 2 ] and thus I k (k — 0, . . . , TV — 2) has 
the physical dimension [L fc ~^~ 2 ]. For instance, l\f_ 2 an d ^V-3 are calculated as 

-AAfI^_ 2 = Af(Af - 1),R + 2Af(Af - 2WK - WwR - 2(/V - ifw^JR 
- 2Af(Af - 1) (wRT + WwM^ + A(Af - 1)^,^ 
-4(Ar-l)«l 1 )Wi[L' 4 ], (2-9) 

-12TV/^ 3 = Af(Af - 1)(TV - 2) (w^'l' + 2«,#K) - 6Af(wR + 2^4) 
- (AT - 1) (Af - 2) [(2N - 3)wW lW W-" + 6^J^_>K; 
+ 6(Af - 2) [w^[ W W 2 + (M- 1)<W^ + Mw^A 
+ 12(A/- - 2)w^ 2 w^ 2 - 6(AT - 1)(A/- - 2) [(r.JJ!!,) 2 ^ + w^w^y 
-12(Af-2) W W lW W lW W 2 [L-% (2.10) 

and consist of the terms which are consistent with the dimensional analysis. 

Ideally, one can obtain a general form of Af- fold SUSY systems if one succeeds in express- 
ing all the Af functions w k (k — 0, . . . , Af — 1), which characterize the systems, in terms 
of a single function, say, u and its derivatives u', u", ... by solving the set of the Af — 1 
constraints I k = (k = 0, . . . ,Af — I). If it is eventually the case, we have a set of Af 
functional u k (k — 0, . . . ,7V — 1) such that 

^=4^]^] (A: = 0,...,7V-1). (2.11) 

One of the most important aspects of (2.11) is that each u k (k — 0, . . . ,7V — 1) has the 
same physical dimension as the one of w k . It in particular means that there would be a 
set of transformations w k — > u k which preserve all the physical dimensions. Conversely, 

if we can find a set of dimension-preserving transformations w k — >■ u k , we may solve the 
set of the constraints more easily. In Sections III-V, we will employ this strategy to see how 
drastically we can reduce the complexity of the constraints. 

To solve the constraints to get (2.11) is in principle possible unless some of the constraints 
automatically imply others since we have Af — 1 constraints for the TV unknown functions. 
However, the task would get drastically harder as the integer Af increases. One of the clues 
to circumvent the situation is in (2.4). All the coefficients of derivative operators (except 
for the highest 2/Vth-order and including the lowest Oth-order) in the l.h.s. of (2.4) are 

quadratic forms of w k while the r.h.s. depends on, in addition to the potentials V ± , the 

TV constants Ck (k — 0, . . . , TV— 1) which are absent in the l.h.s. The latter fact indicates the 



existence of M integral constants of any AA-fold SUSY system which would be functionals of 
w \ (j ' — 0) • • • 1-W ~ 1) an d V ± whose physical dimensions are the same as the ones of C k : 



C k = J k [wM,V] [L" 2 ( fe+1 )] (A; = 0,...,JV-1). (2.12) 

They must emerge from the integration of the set of differential equations 

^-J k [w w ,V] =0 [L-< 2fc+3 >] (k = 0,...,Af-l). (2.13) 

The latter equations are another set of constraints. On the other hand, the set of equalities 
Ik = (A; = 0, . . . ,J\f) are the only constraints which come from the condition for A/*-fold 
SUSY. Hence, the differential equations (2.13) must be such equations that hold whenever 
all the conditions I k — (k — 0, . . . ,J\f) are satisfied. This means that all the quantities 
dJk/dq (k = 0, . . . , J\f — 1) in the l.h.s. of (2.13) would be expressible in terms of /• in 

a way such that the identities /!• — (j — 0, . . . , J\f) apparently imply (2.13). The most 
general form of such a kind would be 

^ = EM A "[L- (2fc+3) ] (k = 0,...,tf-l), (2.14) 

q 3=0 

where L k j are all linear differential operators whose coefficients consist of only w k , V ± , 
and their derivatives. Each L k j must has the physical dimension [L^ _2fe_ - 7-1 ] since the ones 

of dJ k /dq and /■ are [L~( 2fc+3 )] and [lA -2- -?], respectively. 

To see the validity of the above argument, let us consider the constant Cq. From (2.8), 
we immediately know the form of Jq as 



C = J [wW, V] = -V- - jfwW-2 + ^f< V i + j^ W-J [ L_2 1- ( 2 - 15 ) 

Hence, the differential equation which leads to the latter equation is 



n - dJ ° - v-' l w [M] ' + l w [Af] " + l w [M] w W] ' [T~ 3 1 (1 \ti\ 

We then check by using (2.7) that the r.h.s. of (2.16) can be in fact expressed in terms of 

/fas 

d J o _ l .„w m _ l T w\ [T -3i ( o 1 7 ^ 

which indeed has the form of (2.14) with 

L oN = JfvW-i [L -1 ]. L o^-i = -jf [L°], L 0k = 0(k = 0,...,Af-2), (2.18) 

all having the correct physical dimensions [Li "• ? ' -1 ] for L j (j = 0, . . . , A/"). 

In practice, we already solved the two conditions (2.7) to obtain the general form of 
V ± as (2.8), and thus we can totally eliminate V ± in the remaining conditions I k = 

6 



(k — 0, . . . , M — 2). As a result, the remaining integral constants C k (k = 1, . . . , M — 1) 
would be functionals of only vn (j — 0, . . . , M — 1): 

C fe = J fc [«;^,V[Ti;^]]:= J fc [ W M] [L" 2 ( fc+1 )] (fc = 1, . . . ,Af - 1). (2.19) 

Accordingly, d.J k /dq (k — 1, . . . ,Af — 1) would be expressed in terms of the remaining /■ 
(j = 0,...,Af-2)as 

o= dJ, = ^ 2 Lfc ^ [L -( 2 , + 3) ] (fc = i > ... j ^_i ) . (2.20) 

q 3=0 

We will later see in Sections III-V that the above analysis is actually valid and helps us to 
obtain general forms of W-fold SUSY for Af = 2, 3, and 4. 

To summarize, the existence of Af constants C k (k — 0, . . . , Af — 1) in the r.h.s. of (2.4) 
has the direct relation to the existence of M constraints I k = (k = 0, . . . , Af — 1) which 

are differential equations after eliminating the algebraic constraint ijy =0. The physical 
dimensions of the relevant quantities are 

/f 1 [L k -*-% tifl [L^], J k [L- 2 ^ 1 )], L kj \i/ f - n ~'- 1 ]. (2.21) 

B. Equivalent Classes of Linear Differential Operators 

The existence of Af — 1 constraints /]• — (j — 0, . . . , J\f — 2) after the determination 
of the potential pair (2.8) results in the existence of null operators associated with the 
A/"-fold SUSY system under consideration. Let us first introduce a linear space of linear 
differential operators, denoted by /C^, whose coefficients are all functionals of only w k 
(k = 0, . . . , Af - 1). Let K tj [wW] e KW (i = 0, 1, 2, . . .) and define a set of functionals 

ifV^l by 

fW[ w W] = J^K tJ [w^]lf ] . (2.22) 

3=0 

We then define a subspace of K\ \ denoted by K, , which consists of linear differential 
operators whose coefficients are all given by f\ introduced in (2.22). That is, Kq G )C 
means that there exists a set of linear differential operators i^[W '] G /C' ' such that 

Ko = £ if 1 [wM]d t = Y: ( E ^ [- [A/] ]^) <*■ ( 2 - 23 ) 

j j \j=0 / 

It is obvious by definition that any element of K, is a null operator so long as all the Af-fo\d 
SUSY constraints /■ — (j — 0, . . . , Af — 2) are satisfied. Hence, the linear space in which 

an Af- fold SUSY system is considered is actually the quotient space /C' V^o ■ ^ naturally 
leads us to introduce an equivalence class of linear differential oerators in /C' I We shall 
say that two linear differential operators Li,L 2 G /C'^ belong to equivalent class associated 



with Af-fold supersymmetry and express the equivalence as L 2 ~ L\ if L 2 — L\ E JC Q . Any 
equality between operators L 2 = L\ appeared in an A/"-fold SUSY system should be thus 

regarded as an equivalent relation L 2 ~ L\ of the latter equivalent class. In particular, any 
explicit expression for a specific operator such as H ± and Pj^ should be considered as a 
representative of it with respect to the equivalent class. In what follows, we will employ the 

equivalence relation ~ only when we would like to stress that the left and the right hand 
sides of the formula under consideration is identical if and only if (some of) the constraints 
/M = (j = 0, • • • , N - 2) are satisfied. 

III. 2-FOLD SUPERSYMMETRY 

It is quite instructive to see how the general consideration in the previous section make 
sense in the case of 2-fold SUSY though its general form was already obtained by the direct 
integrations of the constraints [2, 11, 12]. Components of 2-fold supercharges are given by 

p- = d 2 + wxd + w , P 2 + = d 2 - w x d + w - w[, (3.1) 

where and hereafter we shall omit the superscript [A/] of w k etc. for the simplicity unless 
the omission would not cause any ambiguity or confusion. The condition for 2-fold SUSY 
P 2 H~ — H + P 2 ~ = is satisified if and only if the following three equalities hold: 

V + - V- = w[, (3.2) 

w" + 2w' + AV~' - 2w 1 (V + - V~) = 0, (3.3) 

w » + 2V-" + 2 Wl V~' - 2w (V + - V~) = 0. (3.4) 

Substituting (3.2) into (3.3) and (3.4), and integrating the resulting equation from (3.3), we 
obtain 

4U+ = 3w[ - 2w + (u>i) 2 - 4C , (3.5) 

AV- = -w[ - 2w + (w,) 2 - 4C , (3.6) 

— 4/ = w'l' — wiw'i — 2(w[) 2 + Aw[wq + 2wiw' — 2(wi) 2 w[ =0, (3.7) 

where Cq is an integral constant. To integrate the third equation (3.7), we shall first make a 
dimension-preserving transformation wq — > uq. We choose it such that it will convert simul- 
taneously both the pairs P 2 and V ± into symmetric forms. The most general transformation 
of polynomial type preserving the physical dimension [L -2 ] of Wq would be 

w = u + -w[ - a (^i) 2 [L~ 2 ], (3.8) 

where Uq [L -2 ] and a is a dimensionless parameter. We note that the latter transformation 
indeed renders both the pair of 2-fold supercharge components and the pair of potentials of 
symmetric forms as 

p± = d 2 T w x d + u Q - a { Wl ) 2 T 7>w[, (3.9) 

AV ± = -2u Q + (2a + l)(wi) 2 ± 2w[ - 4C . (3.10) 



With the transformation (3.8), the condition (3.7) reads as 

-47 = w'C + Aw[u + 2w!u' - 2(4a + l)(w 1 ) 2 w[ = [L" 4 ]. (3.11) 

Hence, it gets simplest when 

a = -1/4. (3.12) 

The next task we should do is to construct the total differential dJi/dq in (2.20). We note 
that dJi/dq and Iq have the physical dimensions [L -5 ] and [L~ 4 ], respectively. Hence, the 
operator L\q defined in (2.20) in this case must have the physical dimension [L -1 ]. Except 
for the differential operator d/dg, there is essentially only one multiplicative operator of 
polynomial type which have that dimension, namely, L\q oc w± [L -1 ]. In fact, we can easily 
check that W\Iq is of a total differential form and thus we put 

dJ 

16— = SwJo = 2w 1 w'{' + Sw^uo + 4(u? 1 ) 2 w' = [L~ 5 ]. (3.13) 

dq 

The latter differential equation is integrated to yield 

I6J1H = 2w i w i - K) 2 + 4(u>i) 2 « = lQC i [ L ~ 4 ]> ( 3 - 14 ) 

where C\ is another integral constant having the correct physical dimension [L~ 4 ] listed in 
(2.5). Hence, we can express uo, and thus wq as well, in terms of w± as 

w[ { Wl ) 2 w'i K) 2 4d , . 

u ^ w °-Y-^--^ + k^ + W- (3 ' 15) 

Substituting it into (3.9) and (3.10), we finally get the general form of 2-fold SUSY systems 

as 

^+ 2 o2 (wi) 2 w'{ K) 2 4Ci w[ 

1 4 2wi 4(wi) 2 (wx) 2 2 v ' 

V ^^ + 4^-8KF"KF ± T- C7 °- (3 ' 17) 

Finally, products of the components of 2-fold supercharges P^P^ 1 are calculated as 

(3.18) 
(3.19) 

where Jo and C\ are given by (3.7) and (3.14), respectively. Hence, we obtain the equality 
(2.4) for M = 2 as an equivalent relation associated with 2-fold SUSY: 

P?P£ ~ ^[( H± + Co) 2 + Ci]. (3.20) 



2 Pi '' 


= 4[(#+ + Co) 2 + Ci] - 2/ , 


2 Pi '' 


= 4[(iJ- + Co) 2 + C 1 ]+2/ , 



IV. 3-FOLD SUPERSYMMETRY 

Next, we shall reexamine the case of 3- fold SUSY, which was once investigated briefly in 
Ref. [13], by utilizing our general analysis. Components of 3- fold supercharges are given by 

p- = d 3 + w 2 d 2 + wid + w , t s 

(4 1) 
p+ = -d 3 + w 2 d 2 - ( Wl - 2w' 2 )d + w - w[ + w'l y ' ; 

The condition for 3-fold SUSY P%H~ — H + P^ = is satisfied if and only if the following 
four equalities hold: 

V + - V- = w' 2 , (4.2) 

w'l + 2w[ + QV - 2w 2 (V + - V~) = 0, (4.3) 

w'i + 2w' + W" + 4w 2 V-' - 2 Wl (V+ - V~) = 0, (4.4) 

w » + 2V'"' + 2w 2 V~" + 2 Wl V - 2w Q {V + - V~) = 0. (4.5) 

Substituting (4.2) into (4.3)-(4.5) and integrating the resulting equation from (4.3), we 
obtain 

6U+ = 5w' 2 - 2wi + (w 2 ) 2 - 6C , (4.6) 

6V~ = -w' 2 -2w 1 + (w 2 ) 2 - 6C , (4.7) 

-6/i = 3w' 2 " + 3w'( - 6w[, - 4w 2 w' 2 ' - 6{w 2 ) 2 + 6w>i + 4w 2 w[ - 4{w 2 ) 2 w 2 = 0, (4.8) 
— 6Jo = w 2 " + 2w"' — 3w'q — w 2 w 2 — 6w' 2 w 2 + w 2 W\ + 2w 2 w'[ + <ow 2 wq + 2w\w' x 

- 2{w 2 ) 2 w 2 - 2w 2 {w 2 ) 2 - 2w 2 w 2 w 1 = 0, (4.9) 

where Cq is an integral constant. To integrate the remaining equations (4.8) and (4.9), we 
shall first construct a set of dimension-preserving transformations Wk — > Uk (k = 0, 1) which 
will convert simultaneously both the pairs P 3 and V ± into symmetric forms, like (3.8) in 
the case of 2-fold SUSY. The most general transformations of polynomial type preserving 
the physical dimensions [L fc ~ 3 ] of Wk would be 



wi = 6-ui + w 2 - ai(w 2 ) 2 [L 2 ], 

wq = uq + 3u[ - f3iw 2 - a\w 2 w' 2 - 6/3 2 w 2 ui - /3 3 (w 2 )' s [L 



( 4 - 10 ) 



where u^ [L fe_3 ] and oto, (5k (k = 1,2,3) are all dimensionless parameters. They actually 
render both the pair of 3-fold supercharge components and the pair of potentials of symmetric 
forms as 

p± = T d 3 + w 2 d 2 t [6«i - a^w^ 2 t w' 2 ]d 

+ u - /3iW 2 - Q/3 2 w 2 u l - f3 3 (w 2 ) 3 =F (3^ - aiW 2 w' 2 ), (4.11a) 

6V ± = - 12ui + (2ai + l)(w 2 ) 2 ± 3w 2 - 6C . (4.11b) 



10 



With the transformations (4.10), the two constraints (4.8) and (4.9) are equivalent to the 
following new set of conditions: 

-3/i = 97i 

= 3O0i + 1)< - 3w + 18(/3 2 + l)w' 2 ui + 6(3/3 2 + 2)w 2 u' l 

- {la 1 - 9/3 3 + 2){w 2 fw' 2 = [L~ 4 ], (4.12) 

3/ = - 9(2/ - dh) 

= 3m'/' - («i - l)u> 2 u> 2 " - 3( 2 l3i + «i + 1)^>2 + 6w' 2 u 
+ 72wi«i - 12(2«i + 3/3 2 + l)w 2 w 2 Mi - 12ai(tu 2 ) M 
+ 2[2( ai ) 2 + ai - 3/3 3 ](u> 2 ) 3 u>2 = [L- 5 ]. (4.13) 

Hence, they would get simplest when 

ai = l, /?i = -l, /3 2 = -l, #, = 1. (4.14) 

With the latter choice of parameter values, the new set of conditions (4.12) and (4.13) reads 



as 



7i=u(, + 2iu 2 «i = 0[L- 4 ], (4.15) 

7 = u'(' + 2w' 2 u + 2Au l u' l - A(w 2 fu l = [IT 5 ]. (4.16) 

The first equation (4.15) enables us to express Mo i n terms of U\ and w 2 as an indefinite 
integral 

u = -2 fdqw 2 u[ [L- 3 ]. (4.17) 

Next, let us construct the total differential d.J\/dq in (2.20). We note that dJi/dq and T& 
(k = 1,2) have the physical dimensions [L~ 5 ] and [L fe ~ 5 ], respectively. Hence, the operators 
Lik defined in (2.20) in this case must have the physical dimension [L~ fc ]. Thus, L10 is just 
a dimensionless constant while L n oc w 2 if we restrict L^ to multiplicative operators of 
polynomial type. Indeed, we can easily find that one of the latter choices leads to a total 
differential 

-4—^ = 7 + 2w 2 h = < + 2(w' 2 u + w 2 u' ) + 2Au 1 u' l = [L" 5 ], (4.18) 

dq 

which can be easily integrated as 

-4JiH = u'[ + 2w 2 u Q + 12( Ml ) 2 = -4Ci [L" 4 ], (4.19) 

where C\ is another integral constant having the correct physical dimension [L~ 4 ] listed in 
(2.5). From (4.16) and (4.19), we can express Uq in terms of u\ and w 2 without recourse to 
any indefinite integral as 

< + 24mi< - A(w 2 ) 2 u[ = < + 12(m 1 ) 2 + 4C 1 
2w' 2 2w 2 

To construct the second integral C 2 of 3-fold SUSY systems, we first note that d J 2 /dq and 7 fc 
(k = 0, 1) have the physical dimensions [L~ 7 ] and [L fe_5 ], respectively. Hence, the operators 
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Z/2fc defined in (2.20) in this case must have the physical dimension [L~ fe ~ 2 ]. Thus, candidates 
for L 2 o are w' 2 , u\, and (W2) 2 , while those for L 2 i are w 2 , u[, u , W2w' 2 , W2U1, and (W2) 3 , if 
we restrict L 2 k to multiplicative operators of polynomial type. With the choice of L 2 o = 
and L 2 i oc -u ) we can construct a total differential as 

dJ 

A—— = u oh = uqu' q + 2w2u[uo 
dq 

= u u' - u\v![ - 12(m 1 ) 2 m / 1 - Adu^ = [L~ 7 ], (4.21) 

where (4.19) has been used. The latter relation is indeed easily integrated as 

8J 2 [w] = (u ) 2 - K) 2 - 8(«i) 3 - 8Ci«i = 8C 2 [L~% (4.22) 

with another integral constant C2 having the correct physical dimension [L~ 6 ] listed in 
(2.5). With the use of (4.22) we can express u solely in terms of U\. Then, substituting the 
obtained expression for uq into (4.15), we can also express W2 solely in terms of u±. Hence, 
we can eventually have an expression for V ± and P 3 in terms of only a single arbitrary 
function u\. However, the latter expression is relatively complicated and thus it would be 
more convenient in practice to express them in terms of two of the three functions W2, u±, 
and uq. If we eliminate uq in (4.11) by using (4.20), we have the expression in terms of W2 
and Mi as 

P£ ~ =F d 3 + w 2 d 2 t [6wi - {w 2 f T w' 2 ] d 

+ w' 2 ' + 6w 2Ul - (w 2 ) 3 ~^-~ ^^- - — =F (3wi - w 2 w' 2 ), (4.23a) 

2w 2 w 2 w 2 

V ± = - 2«i + ^(w 2 ) 2 ± \w' 2 - C . (4.23b) 

On the other hand, if we eliminate W2 in (4.11) by using (4.15), we have the expression in 
terms of u\ and uq as 



**-**-£?* 



6«1 - T7^ ± 



' ^2 / „,// „,"„,' 

, "i"o 



u o) , / u Q U-^Uq 



4K) 2 \2u' x 2(«i 



l\2 



d + u --^- + 



u"> u"u" 



2u[ (u 



l\2 



u>» K) 2 3»A, «) 3 _f 3u , « , <Ko) 2 \ ,j,j , 
U«) 2 K)3 „i J M ° + 8( Wl ) 3 T V x 4K) 2 + 4K)3 J ' ( ' " j '" 



1 '* ;> - 2 " ; + ^C =F (^ - ^4 ) - Co- (4.24b) 



K) 2 + V 4 «i 4 K^ 2 



Finally, products of the components of 3-fold supercharges P 3 T P 3 are calculated as 

P 3 -P 3 + = 8[(77+ + Co) 3 + Ci(77+ + C ) + C 2 ] + 3/i<9 2 + 2[(29 + u> 2 )/i - 7 ]d 

+ (2<9 2 + 2w 2 <9 - 2w 2 + Wl )h - 2(0 + w 2 )I , (4.25) 

P+Pi = 8[(H~ + Co) 3 + C X {H- + C ) + C 2 ] - 3/i9 2 - 2[(d - w 2 )7i + 7 ]9 

-w l I 1 -2(d-w 2 )I , (4.26) 

where 7i, 7 , C±, and C2 are given by (4.8), (4.9), (4.19), and (4.22), respectively. Hence, we 
obtain the equality (2.4) for M = 3 as an equivalent relation associated with 3-fold SUSY: 

P^P± I 8[(77 ± + Co) 3 + C 1 (77 ± + Co) + C 2 ]. (4.27) 
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V. 4-FOLD SUPERSYMMETRY 

In this section, we shall study the case of 4-fold supersymmetry. Components of 4-fold 
supercharges are given by 

P 4 ~ = d A + w 3 d 3 + w 2 d 2 + Wl d + w , 

p+ = d 4 - w 3 d 3 + {w 2 - 3w' 3 )d 2 - (wi - 2w 2 + 3w 3 )d + w - w[ + w 2 - w 3 . ( °' 

The condition for 4-fold supersymmetry P^H~ — H + P^ = is satisfied if and only if the 
following five equalities hold: 

V + - V~ = w 3 , (5.2) 

w'l + 2w 2 + 8V - 2w 3 (V + - V~) = 0, (5.3) 

w' 2 + 2w[ + \2V~" + 6w 3 V~' - 2w 2 (V + - V~) = 0, (5.4) 

w'l + 2w' + 8V-'" + 6w 3 V~" + Aw 2 V~' - 2 Wl (V + - V~) = 0, (5.5) 

< + 2V-"" + 2w 3 V" + 2w 2 V-" + 2 Wl V~' - 2w (V + - V~) = 0. (5.6) 

Substituting (5.2) into (5.3)-(5.6) and integrating the resulting equation from (5.3), we 
obtain 

SV + = 7w' 3 - 2w 2 + (w 3 ) 2 - 8C , (5.7) 

SV~ = -w 3 - 2w 2 + (w 3 ) 2 - 8C , (5.8) 

-8I 2 = 6w 3 + 8w% - 8w[ - 9w 3 w 3 - 12{w 3 ) 2 + 8w' 3 w 2 + Qw 3 w 2 - Q{w 3 ) 2 w 3 = 0, (5.9) 
-8/i = Aw 3 " + 8w 2 - Aw" - 8w' - hw 3 w 3 - 2Aw' 3 w'l + 2w' 3 w 2 + Qw 3 w 2 

+ 8w' 3 wi + Aw 2 w 2 — 6(w 3 ) 2 w 3 — 6w 3 (w 3 ) 2 — Aw 3 w' 3 w 2 = 0, (5.10) 

-8/ = w 3 " + 2w 2 - Aw'v - w 3 w 3 " - 8w 3 w 3 - Q(w 3 ) 2 + w 3 w 2 + 2w 3 w 2 
+ w 3 wi + 8w 3 wo + 2w 2 w 2 + 2w 2 wi — 2{w 3 ) 2 w 3 — 6w 3 w 3 w 3 
— 2w 3 w 3 w 2 — 2{w' 3 ) 2 w 2 — 2w 3 w' 3 wi = 0. (5-H) 

where Co is an integral constant. To integrate the remaining equations (5.9)-(5.11), let us 
first look for a set of dimension-preserving transformations Wk — > Uk (k = 0, 1, 2) which will 
convert simultaneously both the pairs P 4 and V ± into symmetric forms as in the cases of 
2- and 3-fold SUSY. The most general transformations of polynomial type preserving the 
physical dimensions [L fe ~ 4 ] of w^ would be 

w 2 = u 2 + -w' 3 - a ± (w 3 ) 2 [L~ 2 ], (5.12) 

w\ = Mi + u 2 - 6iW 3 - 2aiw 3 w 3 - 2 w 3 u 2 - /3 3 (w 3 ) 3 [L~ 3 ], (5.13) 

/ // I Pi \ III II / / \2 I 2 / \i 

w = u + -u x - 7i^ 2 - I — + - I w 3 - -f 2 w 3 w 3 - 73(^3) - y ( w 3^ 2 ) 

33 
- 74W3M1 - 75(«2) 2 ^(w 3 ) 2 w 3 - 7 6 (w 3 ) 2 M2 + l7(w 3 ) 4 [L~ 4 ], (5.14) 

where u^ [L fe ~ 4 ] and a 0) Pk (k = 1, 2, 3), and 7^ {k — 1, . . . , 7) are all dimensionless param- 
eters. They indeed make both the pair of 4-fold supercharge components and the pair of 
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potentials symmetric as 



Ft = d 4 T w 3 d 3 + 



3 

u 2 -ai(w 3 ) 2 T ^ w 3 



d 2 T 



ui - f3iw'l - /3 2 w 3 u 2 - As(w 3 ) c 



=F (% - 2«iw 3 w 3 ) 



d + u - 7i«2 - 72W3W3 _ 73(^3) - 74W3M1 - 75(^2) 



7 6 (w 3 ) 2 m 2 - 77(^3) 4 =F 



-Ul - 







1 

- Iffl 
2 4 



#2/ w 3/3-j 

:; - — {W3U2) — (w 3 ) W 3 



8\/ d 



2u 2 + (2«! + l)(w 3 ) 2 ± 4u> 3 - 8C - 



(5.15a) 
(5.15b) 



With the transformations (5.14), the remaining three constraints (5.9)— (5.11) are equivalent 
to the following new set of conditions: 



h = 4/ 2 

= - (4ft + 9)u> 3 " + 4u[ - 4(/3 2 + l)w' 3 u 2 - (4/3 2 + 3)w 3 u' 2 

+ (10ai - 12/3 3 + 3)(w 3 ) 2 w 3 = [L- 4 ], 



(5.16) 



4/i = 4(7! - %) 

= - 2(2 7l + 1)< + 4u(, - (4/3 l74 - 4ai + 4 72 + 9 74 + 2)w 3 < 
+ 2(6ai + 2/?! - 2 72 - 4 73 + 3)w> 3 ' - 4( 74 + l)w 3 «i 
- 2(4 75 + l)u 2 u' 2 - 2(2/3 274 - 1a x - 2/3 2 + 2 74 + 4 7e - l)w 3 w' 3 u 2 



2„ ,/ 



(4/3 274 - 2«! + 3 74 + 4 76 )(u! 3 )% 

[4( ai ) 2 - 10a l74 + 12/3 37 4 + 2 ttl - 3 74 - 4/3 3 + 16 l7 ](w 3 ) 3 w' 3 = [L~ 5 ], (5.17) 



7 = - 4(47 - 2dh + d 2 I 2 ) 

= w' 3 '" + Aw' 3 "u 2 - (4/?x + 3)w 3 'u 2 - 4(471 + l)u> 3 w 2 + u> 3 < + 16w 3 m 
+ 4m 2 mx - (6ai + l)(u> 3 ) 2 w 3 " + 4(2a 1 /3 1 + 2a x + A - 47 2 )u> 3 u>> 3 ' 
+ 8(ai - 2 73 )(w 3 ) 3 - 4(2ai + 4 74 + l)w 3 w 3 «i - 167 5 w 3 (w 2 ) 2 
- 4/3 2 u? 3 w 2 w 2 + 4(2a!/3 2 + /5 2 - 4 76 )(«; 3 ) 2 u> 3 w 2 - A(5 3 (w 3 fu' 2 
+ 4(2a 1 /3 3 + (3 3 - A l7 ){w 3 fw' 3 = [L" 6 ]. 

Hence, they would get simplest when 1 

ax = 3/2, A = -9/4, f3 2 = -1, & = 3/2, 7l = -1/2, 72 = 1, 
73 = 11/8, 74 = -1, 75 = -1/4, 76 = 1/2, 77 = -3/8. 



(5.1* 



(5.19) 



1 Another possible choice could be ot\ = 0, /3i = -9/4, /3 2 = -3/4, /3 3 = 1/4, 71 = -1/2, 72 = -1/2, 
73 = —1/8, 74 = — 1, 75 = —1/4, 76 = 0, 77 = 1/16. With the latter choice, it turns out that the third 
integral C3 admits a simpler expression than (Al) but the expression for the second integral C2 contains 
terms which are linear in u\ 1 in contrast with (5.26) and (5.28). 
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With the latter choice of parameter values, the new set of conditions (5.16)-(5.18) reads as 

I 2 = Au' 1 + w 3 u' 2 = Q [L~ 4 ], (5.20) 

h = u' = [L-% (5.21) 

To = W 3 " + Aw' 3 U 2 + 6w' 3 U 2 + Aw' 3 u 2 + W 3 U 2 + I6W3M0 + 4:U' 2 U t 

- 10(w 3 ) 2 w 3 " - 40w 3 w' 3 w' 3 - 10(w 3 ) 3 + 4w' 3 {u 2 ) 2 + 4w 3 u 2 u 2 

- 2A(w 3 fw' 3 u 2 - 6(w 3 ) 3 u' 2 + 30(u> 3 ) V, = [L~ 6 ]. (5.22) 

The first equation (5.20) enables us to express U\ in terms of u 2 and w 3 as an indefinite 
integral 

4ui = - / dgu> 3 w 2 [L -3 ]. (5.23) 

The second equation (5.21) just means that Uq is a constant. On the other hand, C\. 
(k — 0, . . . ,Af — 1) are the only constants which appear in general form of AA-fold SUSY. 
The physical dimension of u is [L~ 4 ] and thus we can put 

u^ = 2C x [L-% (5.24) 

where C\ is an integral constant having the correct physical dimension [L~ 4 ] listed in (2.5). 
We have now obtained the first integral C\ and thus can skip over the step to construct the 
quantity J\. 

Next, to construct the second integral C 2 of 4- fold SUSY systems, we first note that 
dJ 2 /dq and Ik (k = 0,1,2) have the physical dimensions [L~ 7 ] and [L fc_6 ], respectively. 
Hence, the operators L 2 k defined in (2.20) in this case must have the physical dimension 
[L _fc_1 ]. Thus, candidates for L 20 are w 3 , those for L 2 \ are w' 3 , u 2 , and (w 3 ) 2 , while those 
for L 22 are w 3 , u 2 , u±, w 3 w' 3 , w 3 u 2 , and (w 3 ) 3 , if we restrict L 2k to multiplicative operators 
of polynomial type. With the choice of L 20 oc w 3 and L 2 \ = L 22 = 0, we can construct a 
total differential as 

d / 

-128—^ = 2u? 3 / = [L" 7 1, (5.25) 

dq 

where we have omitted the explicit expression. The integration of the latter yields 

-I28J2H = 2w 3 w 3 " - 2w 3 w 3 + (w 3 ) 2 - 16(«i) 2 + 8w 3 w 3 u 2 - A(w' 3 fu 2 

+ 4w 3 w' 3 u' 2 + 2{w 3 fu' 2 ' - 20(u> 3 ) Vj' - 10(w 3 ) 2 (w' 3 ) 2 + 4(w 3 ) 2 (u 2 ) 2 

- I2{w 3 fu 2 + 10(w 3 ) 6 + 32d(w 3 ) 2 = -128C 2 [L- 6 ], (5.26) 

where (5.24) has been applied, and C 2 is another integral constant having the correct physical 
dimension [L -6 ] listed in (2.5). From (5.22) and (5.24), we can express U\ in terms of u 2 and 
w 3 without recourse to any indefinite integral as 

r mil a HI r II I A I II HI 1 1 nl \2 III 

Ui = [—w 3 — Aw 3 u 2 — bw 3 u 2 — 4:W 3 u 2 — w 3 u 2 + W{w 3 ) w 3 

+ A0w 3 w' 3 w 3 + 10(-u; 3 ) 3 — Aw' 3 (u 2 ) 2 — Aw 3 u 2 u 2 + 2A(w 3 ) 2 w' 3 u 2 

+ Q(w 3 ) 3 u' 2 - 30(w 3 ) 4 w> 3 - 32C 1 w' 3 ]/(Au' 2 ). (5.27) 
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Instead, using (5.26), we can express u\ in terms of u 2 and w 3 as a solution to the following 
quadratic equation 

16(«i) 2 = 2w 3 w' 3 m - 2w' 3 w 3 " + {w' 3 ) 2 + 8w 3 w'£u 2 - 4(w 3 ) 2 w 2 + 4w 3 w' 3 u' 2 
+ 2{w 3 ) 2 u' 2 ' - 20(w 3 ) 3 w> 3 ' - W(w 3 ) 2 (w' 3 ) 2 + A(w 3 ) 2 (u 2 ) 2 
- 12(w 3 ) 4 u 2 + W(w 3 ) 6 + 32d(w 3 ) 2 + 128C 2 . (5.28) 

To obtain the third integral C 3 of 4-fold SUSY systems, we note that dJ 3 /dq and Ik (k = 
0,1,2) have the physical dimensions [L~ 9 ] and [L fc ~ 6 ], respectively. Hence, the operators 
L 3 k defined in (2.20) in this case must have the physical dimension [L~ fc ~ 3 ]. With a similar 
dimensonal analysis, we find that the choice of L 30 oc w 3 + Au\ + 2w 3 u 2 — 2(w 3 ) 3 and 
L 31 = L 32 = leads to a total differential: 

dJ 
512-1 = [ W >1 + 4ui + 2w 3 u 2 - 2(w 3 ) 3 ]I = [L~ 9 ], (5.29) 

aq 

where we have omitted the explicit expression. The integral of the latter equation reads as 
J 3 [w] = C 3 with another integral constant C 3 having the correct physical dimension [L~ 8 ] 
listed in (2.5). The explicit formula for J 3 is given by (Al) in Appendix A. By using (at 
least) two of the equalities (5.27), (5.28), and (Al), we can eliminate Ui to obtain the relation 
between w 3 and u 2 . Hence, we are now, in principle, able to express a 4-fold SUSY system 
solely in terms of a single arbitrary function. As in the case of 3-fold SUSY, however, it 
would be more convenient for a practical purpose to have an expression in terms of two of 
the four functions w 3 , u 2 , ui, and u . For example, if we eliminate u\ and u in the system 
by using (5.27) and (5.24), respectively, the system (5.15) with the parameter values (5.19) 
can be represented in terms of w 3 and u 2 as (A2). 

Finally, products of the components of 4-fold supercharges PfP^ are calculated as 

4 
P 4 P+ = 16[(H + + C ) 4 + C 1 (H + + C ) 2 + C 2 {H+ + C ) + C 3 ] + ]T f+[w]d\ (5.30) 

4 
P+Pi = 16[(H- + Co) 4 + C^H- + C ) 2 + C 2 (H~ + C ) + C 3 ] + ^ fr[w}d\ (5.31) 



j=0 



where Ci, C 2 , and C 3 are given by (5.24), (5.26), and (Al), respectively. Both of the second 
terms in the above are elements of /C introduced in Section II B and have the form of 
(2.23). The explicit forms of /^ are given by (A3) and (A4). Hence, we obtain the equality 
(2.4) for M = 4 as an equivalent relation associated with 4-fold SUSY: 

P^P± I 16 [(#* + Co) 4 + C 1 (/f ± + Co) 2 + C 2 (iJ ± + Co) + C 3 ]. (5.32) 

VI. DISCUSSION AND SUMMARY 

In this work, we have clarified general structure of J\f-fold SUSY systems by considering 
dimensional analysis and introducing the equivalent classes of linear differential operators 
associated with them. We have then shown that the latter general consideration is in fact 
effective in constructing the most general AA-fold SUSY systems and their integral constants 
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for M = 2, 3, and 4. Application to systems for H > 4 would be straightforward and the 
problems would get more transparent even though still remain highly complicated. Finally, 
some remarks on the future issues are in order. 

1. Although we have only considered ordinary one-dimensional Schrodinger operators, 
generalization to other operators would be possible. In physical applications, one of the 
interesting extensions is to a quantum system with position-dependent mass for which 
A^-fold SUSY was successfully formulated in Ref. [16]. In the latter case, there is an 
additional freedom of mass function and it is particularly interesting to see how its existence 
would force us to modify or generalize the general considerations made in Section II. 

2. In this work, we have restricted the dimension-preserving transformations (3.8), (4.10), 
and (5.14) to polynomial type, namely, transformations which are polynomials in w k 

(m = 0, 1, 2, . . .). On the other hand, the obtained results such as (3.15), (4.20), and (5.27) 
indicate that the relations among w k are in general expressed by rational functions. 

Hence, we may be able to reduce further the complexity of the conditions for J\f-fold SUSY 
by extending the type of transformations from polynomial to rational function. However, the 
number of admissible forms of rational functions which preserve physical dimension would 
drastically increase. For instance, any sum of rational functions of the form /„ [w]/<7n [w] 
(neZ) 






"f-f+Ei^n'i 



where fn [w] and g n [w] are polynomials in w k having the physical dimensions 

jL.n+fc-A/j and [L n ], respectively, can serve as a (part of) transformation which preserves 

the physical dimension [L fe ~ A/ ] of w k . As a result, we may need additional guidelines to 
restrict the forms of transformations to make an efficient analysis. We are curious to know 
how to get such guidelines systematically for the purpose. 

3. As was pointed out in Section II, an A/"- fold SUSY system is in general only weakly quasi- 
solvable but is not quasi-solvable in the strong sense and thus does not necessarily admit 
analytic local solutions in closed form. In the case of M = 2, it was proved in Ref. [17] that 
type A 2-fold SUSY is a necessary and sufficient condition for a one-dimensional quantum 
mechanical system to have quasi-solvability in the strong sense with two independent analytic 
local solutions. So, it is interesting to see what kind of condition is necessary and sufficient 
for a quantum system to admit three or four independent analytic local solutions. We will 
report on the latter subjects in subsequent publications. 



Appendix A: Formulas 

In this Appendix, we present the lengthy formulas appeared in 4-fold SUSY in Section V. 
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The third integral: 



1024 J 3 [w] = 2w%w' 3 " ~ W) + 8w 3 "«i + 4w 3 w' 3 "u 2 - 4w' 3 w 3 u 2 + 6{w 



H\2„ 



U 2 



— 2w 3 w 3 u 2 + 6w' 3 w 3 u 2 + 2w 3 w 3 u 2 + 32w 3 u 2 u 1 + 24w' 3 u 2 Ui + 8w 3 u 2 Ui 

+ 32u 2 (u 1 ) 2 - 4(w 3 ) 3 w' 3 m + \2(w 3 ) 2 w' 3 w' 3 " - 16(w 3 ) 2 (w' 3 r ) 2 - 24w 3 (w' 3 ) 2 w 3 ' 
+ (w' 3 ) 4 - 80(w 3 ) 2 w^ Ul - 80w 3 (w' 3 ) 2 Ul + lQw 3 w' 3 \u 2 ) 2 - 4(w' 3 ) 2 (u 2 ) 2 
+ 8w 3 w 3 u 2 u 2 + 4(w 3 ) 2 u 2 u 2 - (w 3 ) 2 (u 2 ) 2 + 32w 3 (u 2 ) 2 u 1 - 56(w 3 ) 3 w 3 u 2 

- 20{w 3 ) 2 {w' 3 ) 2 u 2 - 4(u> 3 ) V 2 ' - 64{w 3 ) 3 u 2Ul + 8{w 3 ) 2 {u 2 ) 3 + 40{w 3 ) 5 w' 3 ' 

+ 10(«; 3 ) 4 (u> 3 ) 2 + 48K) 5 Ml - 28{w 3 )\u 2 ) 2 + 3Q{w 3 fu 2 - 15(u> 3 ) 8 + 32C 1 (w 3 ) 2 

+ 2hQC 1 w 3 u 1 + 64C 1 (u> 3 ) 2 w 2 - 32d(w 3 ) 4 + 256(d) 2 = 1024C 3 [L~ 8 ]. (Al) 



P4 and V in terms of w 3 and w 2 : 



P±^ 



d 4 T w 3 9 3 + 



«2~ 2^3) T^3 



<9 2 T 



«'., 



4 W3 -4^ 



w 3 'u 2 

uL 



w' 3 u 2 ' 
uL 



w 3 u 2 
4m' 



5(u> 3 ) 2 -u; 3 " 10u> 3 u> 3 w 3 5(w 3 ) 3 w 3 (-u 2 ) 2 : 6(w 3 ) 2 w' 3 u 2 15(w 3 ) 4 w' 3 



2w 2 
8Ciw' 3 



u,-. 



2u' 



Un 



'/o 



=F («' 2 - 3w 3 w 3 ) 



1 5 

8 + 2 U ' 2 ~ 2 WsW ' 3 ' 



u 2 
f(^) 2 + \(u 2 f 



2u'o 



\{w 3 ) 2 u 2 



+ 1 4m 4 



w 3 w 3 
4u' 



W 3 W 3 U 2 

U 2 
,/ \3 



W3W 3 U 2 

Un 



3„„w 



MX' 5( W3 ) 3 ^ 3 ' 
4 M ; 



2m' 



10(iy 3 ) w' 3 w 3 5w 3 (w 3 ) w 3 w' 3 {u 2 ) 6(w 3 ) w' 3 u 2 15 (w 3 )w' 3 

1 



u 2 
lCiw 3 w' 3 

Un 



2u'o 



a:-, 



+ 2dT 



7 



u 2 
9 



2m' 



-u> 3 + -u> 3 w 2 + -w 3 m 2 - 7(^3) w 3 



^ ±= - |«a + |(wa) 2 ± |ws 



C . 



(A2a) 
(A2b) 



The functions ft[ w ] in (5.30) and (5.31): 



ft 
ft 
ft 

ft 



-4/ 2 , /+ = 4/ 1 -4(39 + W3 )/ 2 , 

- 4/ + {88 + 3iu 3 )/i - i[289 2 + 15«; 3 <9 - 9w' 3 + 4w 2 + {w 3 ) 2 ]I 2 , 

- 2(29 + w 3 )I + 2(3<9 2 + 2w 3 d - 2w 3 + w 2 )h 

- 2[4<9 3 + 3w 3 d 2 - 2(2w 3 - w 2 )d - 4w 3 + w 1 - 2w 3 w 3 ]I 2 , 

- 2(d 2 + w 3 d - 2w' 3 + w 2 )I + [2d 3 + 2u? 3 <9 2 - 2(2w' 3 - w 2 )d - 4w' 3 ' 

- 2w 3 w']h [32<9 4 + 32w 3 d 3 - 32(2w' - w 2 )d 2 - 4(3lw' 3 - 2w 2 
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w 1 
6wi 



+ 14u? 3 u? 3 )<9 - 60w 3 + 8w 2 + 16w - 68w 3 w' 3 + 7{w' 3 ) 2 - 4w' 3 w 2 + 8w 3 w 2 
- 4(w 2 ) 2 - 22(w 3 ) 2 w' 3 + 4(w 3 ) 2 w 2 - (w 3 ) 4 ]I 2 , 



(A3) 
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and 



/- = 4/ + (Ad - 3w 3 )h + l -[Ad 2 - 3w 3 d - 3w' 3 + Aw 2 - (w 3 ) 2 ]I 2 , 



/ 4 - = 4/ 2 , f~ = Ah + A{d - w 3 )I 2 , 

/f - 2(29 - w 3 )/ + 2(<9 2 - w 3 <9 - u> 3 + w 2 )h + 2(w'l + 2w' 2 - Wl - 2w 3 w' 3 )I 2 , 

f = 2(8 2 - w 3 -w' 3 + w 2 )I + [w'i + 2w' 2 - Wl - 2w 3 w' 3 )h - —[A{w'i + 2w' 2 

lb 

- 2iwi - 2w 3 w 3 )9 - 12w 3 - 2Aw 2 + 16w + 28w 3 w 3 + 23(w 3 ) 2 - Aw' 3 w 2 
+ 8w 3 w' 2 - A(w 2 ) 2 - 6(w 3 ) 2 w 3 + A(w 3 ) 2 w 2 - (w 3 ) A ]I 2 , 

where I 2 , I±, and I are given by (5.9), (5.10), and (5.11), respectively. 
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